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Review on string breaking 

- the query in quest of the evidence - 
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Considerable progress has been achieved recently in the observation of string breaking within non-Abelian Higgs 
models, by use of multi-channel methods allowing for broken string states. Similarly, in pure gauge theory this 
approach has been shown to reveal string breaking for color charges in the adjoint represenation. For QCD with 
dynamical fermions, one needs substantial progress in noise reduction, however, to render such techniques viable. 



1. INTRODUCTION 

Twenty years after the demonstration of con- 
finement in quenched QCD simulations (!]] we 
have beautiful evidence of the formation of colour 
flux tubes between static colour charges ||. As 
to the verification of their fission, it has been an- 
ticipated ever since that full QCD calculations 
should reveal string breaking (SB), in form of a 
screening behaviour in the static potential V(R). 

While studies on Polyakov loop correlators in 
the temperature range .68 < T/T c < .94 § 
reviewed at LATTICE98 {§ did indeed provide 
clear evidence for such phenomenon (c/. also 
the arguments presented in ||), the head-on ap- 
proach for detection of SB signals from mere Wil- 
son loop studies at T = has not proven success- 
ful ||. This situation is illustrated in Fig. |l| e.g., 
with updated T%L data 0. 

But why bother about SB if one expects it any- 
how? The motivations to focus on the challenge 
are obvious: (a) the i?-asymptotia of V(R) is the 
most obvious infrared quantity that we can tar- 
get on the lattice and hence important to be un- 
der control, (b) SB can teach us both about the 
techniques to handle mixing problems and hadron 
decays like T — > BB on the lattice ||. Lastly, 
understanding hadronization in the static inter- 
action scenario of full QCD will certainly help 
to further substantiate previous quenched results 



on confinement |i0|| . For we should remember 
that transfer matrix studies yield at best numer- 
ical estimates of upper bounds to V as long as 
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Figure 1. Update of full QCD potential from T%L 
(at (3 — 5.6) |0, in comparison to their quenched 
result at (3 q = 6.0, with Sommer radius r ||. 



we are analyzing asymptotic (in Euclidean time 
T) behaviour of Wilson loops W in terms of 
dlnW{R,T)/dT -> -V(R), at finite T. 
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Torelon masses, k = 0.1575 
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Figure 2. Plateau formation in V^//(2.4ro) at = 
5.6 and n sea — .1575, from T%L configs. 




Figure 3. Effective torrelon masses from SESAM 
configurations at n sea — .1575 E3. 



2. VISIBILITY: 5 LATTICE SPACINGS 

Quenched analyses of static potentials rely - 
apart from high statistics - on the combined ac- 
tion of three signal enhancing techniques : (a) 
smearing the spatial links in order to suppress ex- 
cited state contaminations, (b) analytic multihit 
noise reduction on the time- like links, (c) exploit- 
ing loop averaging effects on each configuration 
by shifting and measuring W all over the entire 
lattice. 

In full QCD we are in lack of high statistics 
samples of vacuum configurations (for cost rea- 
sons) and moreover, we have to abandon multi- 
hitting as it becomes unfeasable in presence of 
long range quark loop effects. In Fig.^ I demon- 
strate the quality of the resulting plateau in the 
effective potential, V e //(2.4r , T). Evidently one 
runs into rapidly increasing errors when going be- 
yond T = 5 — 6! 

The lesson to be learnt from Fig.|| is obviously 
that - in order to uncover SB - one better aims 
at achieving an improved overlap with the true 
ground state O], in terms of an earlier onset of 
the plateau, in a region where the resolution of 
the effective potential still suffices, at time sepa- 



rations 3,4, and 5. 

Let me just mention in passing an alterna- 
tive strategy to establish SB: (a) determine in 
full QCD the masses of torelons, i.e. the spatial 
wrappers (of length L) around the lattice 
(b) compute an effective string tension, K, like 
in quenched QCD from a fit to the potential at 
r < 2tq where string breaking is not yet expected 
to set in; (c) torelon effective masses undershoot- 
ing the string energy would provide evidence for 
SB. As illustrated in Fig. | the SESAM data on 
the 16 3 x 32 lattice at their smallest quark mass 
exhibits some weak evidence for SB, by crossing 
the level aL = 16K, though at rather high noise 
level. 

Another aspect which has been brought for- 
ward recently [|l3| is to exploit coarse grained spa- 
tial lattices by help of improved actions. 

3. MULTICHANNEL APPROACH 

The obvious strategy for improving the ground 
state overlap is to enlarge the operator space to a 
multichannel approach comprising hadronization. 
This would be very much in the spirit of multi- 
operator variational ansatze as devised and used 
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largely in the context of glueball studies. In 
the present context strong coupling methods have 
been used to demonstrate that inclusion of 
the QQ qq channel in addition to the static quark- 
antiquark QQ (with fluxtube in between) looks 
like a promising path for uncovering SB. It implies 
however the use of large loops with color partly 
transmitted by insertion of light quark propaga- 
tors (instead of by gauge links only) . This feature 
precludes us from applying easy loop averaging 
over the entire lattice, since light quark propaga- 
tors are too expensive to compute for all source 
points on the lattice. 

4. SB WITHOUT FERMIONS 

Given the problems to deal with fermions refs. 



16 I?]] followed an old suggestion of Bock et al 
18 1 to study the mechanism of SB in an easier 
setting, namely by resorting to non-Abelian Higgs 
fields. 

In the two-channel scenario (see Fig. |J) one 
is faced with a generalized eigenvalue problem of 
the transfer matrix, for evolution from T = to 
T , with eigenvalues X l {T,T Q ) p| 



C{T) Si = Xi{T,T ) C{T )c u 



(1) 



the index labeling ground (i — 0) and excited 
state (i = 1). In a contribution to this conference 
Knechtli et al 123] have presented a high statis- 
tics scaling study of an SU(2) Higgs model in 3+1 
dimensions, with matter field $ in the fundamen- 
tal representation, in the confining phase. As can 
be seen from Fig. |5| they can neatly trace both 
ground and excited states in a scaling regime, 
with a clear-cut gap in between. 
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Figure 4. The two channel transfer matrix C(T) 
in the non-Abelian Higgs model, with due tribute 
to Pisa: the bullets stand for scalar mass inser- 
tions. 



l 
o 
-l 

-2 

-3 - 



j ■ 

excited potential m 



o (3 = 2.4 
* (3 = 2.2 







1 2 
r/r 



Figure 5. Scaling of ground and excited states 
potential in the SU(2) Higgs model ptJ. 
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Figure 6. The overlap of the Wilson loop opera- 
tor, uiq, with the broken string ground state p0|. 



Moreover they show that the overlap function 
of the Wilson loop operator to the true ground 
state, lo^(R), computed according to 

n(T) = (S , C(T) So) = exp (—V n T), (2) 

n 

with O(0) = 1 exhibits a steeply falling in- 
dependence, as shown in Fig. [(]. They illus- 
trate moreover that a naive one-channel analysis 
is manifestly misleading, as it yields an overlap 
estimate, L J^' nmve ; that fakes a good projection 
of the large Wilson loops to the ground state!! 

As yet another check for the efficiency of the 
two-channel approach in revealing SB, two de- 
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Figure 7. SB in adjoint potential compared to 
confinement in static potential in fundamental 
representation in 2+ld SU(2) g.t., from ref. p2|. 
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Figure 8. SB observed for three different excita- 
tion levels of the adjoint string, from ref. [El). 



tailed studies [^T|,^2| were addressed recently to 
the adjoint flux tube formation between static 
color charges in the adjoint representation of 
SU(2) gauge theory in 2+1 dimensions. Such 
strings are not protected from breaking up, even 
in absence of fermions |23|j . In the correlation 
matrix C (see Fig ||) the mass insertions are 
now to be replaced by gluelump operators^, 
gluelump^ = Tr (clover^ a a ) (with color a and 
site x) to be linked by the color flux operator 
Tr(a a U x -, y G h Ul_+y). Fig. | highlights beauti- 
fully the efficiency of the multichannel ansatz in 
revealing, in addition to the ground state, another 
three excited states, along with resolution of their 
level crossings! 

Having convinced ourselves that an opera- 
tor set extended to include physically motivated 
channel mixing is the right way to go, let us return 
to our original problem, i.e. QCD with fermions. 



5. UNCOVER SB IN QCD? 

If we wish to adopt the multichannel strategy 
to the QCD situation, we have to deal with loops 
containing insertion of light quark propagators, 
like in the transition correlator (2|C|1) depicted 
in Fig. 0. 

Fig. [To] [^5| illustrates the noise situation with 
the example of local masses from the transition 




static antiquark static quark 

Figure 9. Transition matrix element from state 
QQ to state QQ qq , with source/sink smearing. 



correlator, (2|C|1). While source smearing helps 
a lot, doing without loop averaging makes things 
hopeless, even at a small value of T, such as 
T = 2(!). This can best be seen by comparing to 
the noise level in the corresponding Wilson loop, 
( 1 1 C 1 1 } ! One should be aware that the noise sit- 
uation is considerably aggravated for (2|C|2), at 
T = 3 or 4! Similar observations have been re- 
ported by Lacock for the MILC collaboration (2(J . 

In order to get noise under control one obvi- 
ously needs the computation of light quark prop- 
agators P(x —►!/) = on source locations 
sampled over a large region, R, of the lattice vol- 
ume. A promising program in this direction has 
been launched recently by Michael et al. |^7j Their 
proposal is to start out from a Gaussian estimate 
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Figure 10. Noise on local masses vs R,&tT= 2: 
small (large) error points refer to Wilson loops 
(transition matrix elements, (2|C|1), cf. Fig. ^|); 
horizontal lines mark the error band for 2mps- 
Analysis done on T%L configs, at n sea — .1575 
and r w 5.9a ||. 



for inverse Dirac operator in terms of a random 
scalar field $. 

M-y 1 = \ j D ® <f;<f,cxp{-$*Af$/2} . (3) 

The estimator introduces additional noise whose 
variance can be minimized by a multihit (analyt- 
ical integration) on <£> over the interior region of 
R, I = R — dR. In praxi, this requires computa- 
tion of an inverse block matrix M _1 (by iterative 
solver) after each update of the <3>-field on the 
boundary dR, as borne out in the replacement: 

®x->v x = / D$ $ z exp{-$*M$/2} 

= -M- X ]M X , Z <£ Z . (4) 

Here M (M) is a block matrix with support C 
I = R -dR(R) and z e dR . 

So far this noise reduction technique has been 
applied successfully to the determination of forces 
between two static-light mesons p8fl . It remains 
to be seen whether it can bring decisive progress 
in the treatment of QCD string fission. 



6. PISA EST OMEN? 

While in confining models without fermions, 
SB has been clearly seen at work, we have at 
this stage no real compelling evidence for string 
breaking from simulation of QCD with dynamical 
fermions. 

In any case, from meditating on Fig. |ll| which 
is a zoom to this years conference poster a su- 
perstitious person might suppose that at present 
Pisa would not offer favourable auspices for the 
occurence of string breaking. So for the time be- 




Figure 11. The Pisa model for b-confmement. 



ing, let me finish by quoting a famous last century 
tourist who felt frustration after climbing (the 
tower?): "Spitze Steine - Miide Beine - Aussicht 
keine - Heinrich Heine". What about Heine be- 
ing wrong? 
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